Ferroelectricity in the Magnetic E-Phase of Orthorhombic Perovskites 
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We show that the symmetry of the spin zigzag chain E phase of the orthorhombic perovskite 
manganites and nickelates allows for the existence of a finite ferroelectric polarization. The proposed 
microscopic mechanism is independent of spin-orbit coupling. We predict that the polarization 
induced by the E-type magnetic order can potentially be enhanced by up to two orders of magnitude 
with respect to that in the spiral magnetic phases of TbMnOa and similar multiferroic compounds. 
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Introduction. The switching of the electric polariza- 
tion P by a magnetic field of a few Tesla discovered in 
TbMn0 3 and TbMn 2 5 2] has ignited enormous in- 
terest in a class of materials that can be termed improper 
magnetic ferroelectrics (IMF's). While there is no intrin- 
sic ferroelectric (FE) instability in the IMF's, P emerges 
due to its coupling to the primary magnetic order param- 
eter. Hence, the FE phase transition coincides with the 
corresponding magnetic transition, and P is very sensi- 
tive to magnetic field induced changes of the magnetic 
state USlSSSSSSIil s y mmetr y imposes rather 
strict conditions on the possible magnetic order parame- 
ters - the magnetic structure must have low enough sym- 
metry in order for the system to form a polar axis. As 
a consequence, the IMF's often have complicated non- 
collinear structures, including spiral and incommensu- 
rate SSSSEQE1E3, while the IMF phases 
with collinear magnetism are rare 0, . Noncollinear 
magnetic structures are stabilized due to either compet- 
ing interactions (frustration) or anisotropies generated 
by spin-orbit coupling, which leads to reduced transition 
temperatures and weaker order parameters. The mag- 
nitude of P (w O.l^C/cm 2 ) is also affected by its weak 
coupling to magnetism. In turn, collinear IMF's may 
prove more promising for future applications as they are 
less prone to the obstacles mentioned above. 

In our present study we turn our attention to the 
collinear E-type magnetic phase that has been observed 
in perovskite manganites 0, 0] and nickelates 0, 0, 
EFl EH ]. First, we show that this is a previously over- 
looked example of an IMF. Second, in contrast to the 
most extensively theoretically studied case of spiral mag- 
netism [22J, |23l |2J], the mechanism responsible for fer- 
roelectricity in the magnetically ordered phase does not 
rely on the presence of anisotropic Dzyaloshinskii-Moriya 
interaction. In our case, P appears due to a gain in the 
band energy of the e g electrons in the double-exchange 
model. We estimate that P can be potentially larger than 
that in the spiral magnets by up to two orders of mag- 
nitude, thus paving a way for IMF's to reach the values 
of P common for other multiferroic systems with better 



FIG. 1: (Color online) (a) Magnetic unit cells of the two E- 
phase domains in HoMn03 corresponding to the (Ei, 0,—P a ) 
and (0, £2, Pa) solutions of ©• The arrows on the Mn atoms 
denote the directions of their spins. The FE displacements are 
not shown, but the direction of P is indicated, (b) The simple 
sine-wave ma gnet ic structure of HoMn03 for Tl < T < Tjv = 
42.2 — 47.5 K (l6|. (c) The E-phase magnetic structure below 
T L . 



ferroelectric properties [251 12i 

In the perovskite manganite family i?Mn03 (space 
group Pbnm), the E-phase was first reported for i?=Ho 
as a result of mag netic structure refinement by neu- 
tron diffraction [16|. The magnetic unit cell of the E- 
phase, forming at low temperature below Tl — 26 — 
29.6 K [III23, is shown in Fig.QIa). The Mn atoms with 
parallel spins form zigzag chains in the a6-plane, with the 
chain link equal to the nearest-neighbor Mn-Mn distance. 
The neighboring zigzag chains in the 6-direction have 
antiparallel spins. Figures ^b) and (c) show that the 
equal-spin E-phase structure can be obtained from the 
simple sine-wave structure by locking-in its modulation 
vector kb = ~ and fine-tuning its phase. The a6-plancs 

are stacked antiferromagnetically (-1 1 — ) along the c- 

direction. Interestingly, Lorenz et al. [271 ] reported on a 
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large magnetic field dependence of the dielectric constant 
below Tl, which may be an indication of the multiferroic 
state. 

Landau theory. Independently of the microscopic 
mechanism, the possibility of a FE state in magnets can 
be examined by considerin g the sy mmetry allowed terms 
in the Landau potential |j TlH El El Ei^ . We define the 
symmetric coordinates corresponding to the E-phase as 

Ei = Si + S2 — S3 — S4 — S5 — Sg + S7 + as, 
E 2 = S 1 -S 2 -S3 + S4-S 5 + S 6 + S 7 -S 8 , (1) 

where Si is the spin of the ith Mn atom in the magnetic 
unit cell, as shown in Fig^a). Since the Mn spins in 
HoMn03 point along the &-axis |l6]. below we consider 
only the ^-components of E12 denoted by -Ei,2- However, 
the expression for the Landau potential derived below is 
valid for any component of Ei,2- E\ and E 2 span an 
irreducible representation of the space group Pbnm cor- 
responding to k = (OiO). The properties of this repre- 
sentation are summarized in Tabled Taking into account 
that P transforms as a polar vector, we obtain the fol- 
lowing form of the Landau potential corresponding to the 
E-phase, 

F = a(Ef + El) + bi(E\ + Elf + b 2 E 2 x El 

+c{E 2 - Ei)P a + d(E 2 - El)E x E 2 P h + -^P 2 , (2) 

where x is the dielectric susceptibility of the paraelectric 
phase, and the other coefficients are phenomenological 
parameters. Minimizing F with respect to P, we ob- 
tain P a = -c X {E\ - El), P b = -d X {El - E 2 )E X E 2 , 
and P c = 0. Therefore, each of the four domains of the 
E-phase [(±£' 1 ,0) and (0,±E 2 )] is IMF with the polar- 
ization along the a-axis and different signs of P a for E\ 
and E 2 . We also notice that the 6-axis component of P 
can be locally present within the domain walls where E\ 
and E 2 coexist [3(J- 

The E-type phase in the nickelates also consists of 
spin zigzag chains that have a different direction in 
the afe-plane and a different stacking along the c-axis 
(+ + ) HEmEl. The corresponding modu- 
lation vector is k = (^0^). The Landau theory analysis 
leads to results similar to the manganite case. For nick- 
elates, the symmetric coordinates E\ and E 2 (not given 
here) are written in terms of 16 Ni spins of the magnetic 
unit cell which is 4 times the crystallographic unit cell. 
It leads to a Landau potential similar to J5J) with P a re- 
placed by Pb in the fourth term and d = 0. Thus, P in 
the E-phase of nickelates is parallel to the 6-axis. 

Microscopic model. To understand microscopically the 
possible mechanism of ferroelectricity in the E-phase we 
use Monte Carlo (MC) simulations to study the ground 
state properties of the following Hamiltonian for mangan- 
ites based on the orbitally degenerate double-exchange 



TABLE I: Matrices of the generators of space group Pbnm in 
the irreducible representation spanned by E\, E2- The space 
group elements are denoted (r\hkl), where r is the identity op- 
eration 1, two-fold rotation 2 a , c , inversion /, or time reversal 
1' followed by the translation t = ha + kh + Ic. 
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model |3~ill33.l33| with one e g -electron per Mn 3+ ion, 

H = — ^ Ci,i+ai^gdi Q di+a/3 + JAF ^ Si • Si +a (3) 
iaa/3 ia 

+A^(QliPi + QftTxi + Q 3 iT z i) + - ^2 KmQmii 
i im 

where d\ a is the creation operator for the e g electron 
on orbital a — x 2 — y 2 (a,),3z 2 — r 2 (b); a = x, y is the 
direction of the link connecting the two nearest neigh- 
bor Mn sites; and Si is the classical unit spin, given by 
the polar angle Qi and azimuthal angle <fii, representing 
the electrons occupying the t 2g orbitals on the i-th Mn 
site. The t 2g spins interact directly through the antifer- 
romagnetic superexchange Jaf>0. Qj = cos If cos + 
sin y sin ^e~ l< ^ i ~'^ i " 1 is the double- exchange factor aris- 
ing due to the large Hund's coupling that projects out the 
e g electrons with spin antiparallel to S.;. Q m \ represent 
the classical adiabatic phonon modes, with stiffnesses K m , 
due to the displacements of the ligand oxygen ions sur- 
rounding the i-th Mn site. The phonon stiffnesses were 
chosen as follows - for the Jahn- Teller modes: k,\ = 2.0, 
«2 = ft 3 = 1-0; for the FE mode kfe = 8.0 and for 
the rest of the modes, k = 10.0. Except for the Jahn- 
Teller modes |U EH , this choice is driven only by the 
efficiency of the MC simulations, and it does not reflect 
the actual frequencies of the phonon modes, which are 
presently unknown. We illustrate below how P can be 
obtained in physical units. The third term in Eq. Q 
is the Jahn- Teller coupling with constant A and the e g - 
orbital operators p ; = d\ a d- la + d\ b dn, r xi = d\ a di b + d\ b d ia 
and r zi = d\ a d ia - d\ b d ib . 

Our improvements over previous approaches to this 
model are the following. First, we explicitly consider 
the dependence of the hoping parameters t^g on the 
angle of the Mn-O-Mn bond tp- la . Taking into account 
only the largest Mn-0 cr-bond contribution, we find [34| 

= *aa = t COS 3 y?, <£ b = i£ b = -<COS^/3, t^ b = 

*ba = -<b = = ^ cos 2 <p/y/3, where t = l(pda) 2 is 
taken as the energy unit hereafter. We neglect the <p de- 
pendence of Jaf since it has a smaller energy scale than 
t. Second, Qmi's are defined such that the elastic energy 
term in J2J is minimal for ip\ & = tpo < 180°. This allows 
us to model the initial structural buckling (GdFe03-type) 
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FIG. 2: (Color online) (a) The starting configuration of a 
Mn-O-Mn bond. Numbers 1-4 enumerate the O atoms sur- 
rounding one Mn. (b) A MC snapshot of the IMF E-phase 
at T=0.01. The ferromagnetic zigzag chain links are shown 
as solid lines. The displacements of the oxygen atoms are 
exaggerated, (c) Left. The local arrangement of the Mn-O 
bonds with disordered Mn spins (full circles). Right. Oxy- 
gen displacements (arrows) within the chains of opposite Mn 
spins (open and crossed circles) in the E-phase (see also 
Refs-laEl). (d) MC results for the polarization at T=0.01 
for different values of Jaf. 

distortion present in the orthorhombic perovskites [3S| . 
In particular, the buckling mode is defined in two dimen- 
sions as Qbuckica = yn - V-a - x i3 + x u - (-l)^ +, «4u , 
where Xik , yik are the displacements of the oxygen atoms 
from their ideal position in the 180° Mn-O-Mn bond, 
u = a cot ^p, and a is the Mn-0 distance in the 180° 
Mn-O-Mn bond [see Fig. [21 a)]. Since the buckling distor- 
tion in the considered perovskites is present at tempera- 
tures well above the magnetic transitions, we consider ipo 
as a fixed parameter of the model. For HoMn03, (po is 
close to 144° [3{|. As shown below, this distortion plays a 
crucial role in generating ferroelectricity in the E-phase. 
Introducing ipo < 180° effectively reduces the symmetry 
of the Hamiltonian @ , although it is still invariant with 
respect to the inversion symmetry centers located at the 
every Mn site. The FE polarization emerges only due 
to the spontaneous symmetry breaking caused by the E- 
phase magnetic order. 

Monte Carlo results. Many results of the model © 
regarding magnetic and orbital order, without the modi- 
fications mentioned above, are reported in Ref. [3^ |. Our 
present treatment essentially confirms those results. The 
ground state phase diagram in two dimensions contains 
the ferromagnetic phase, E-phase, and the Neel-like G- 
phase. Here we focus our attention on the new results 
directly related to the ferroelectricity of the E-phase. A 
typical low-temperature E-phase MC snapshot is shown 
in Fig.[5{b). In accordance with the absence of spin-orbit 
interaction in Hamiltonian ©, it is invariant with re- 
spect to collective spin rotations, therefore the preferred 
spin direction is chosen randomly in our MC simulations. 
However, the ferromagnetic zigzag chains are clearly es- 



tablished in the ground state. As is seen from the snap- 
shot, the double-exchange physics plays a crucial role in 
the formation of the FE state. Due to the factor Ci i +a 
the electron hopping is prohibited between Mn atoms 
with opposite ti g spins. Hence, the displacement of the 
corresponding oxygen atom perpendicular to the Mn-Mn 
bond (these displacements are not Jahn- Teller active) is 
only due to the elastic energy, which favors the bond an- 
gle ipo. On the contrary, homing is allowed along the 
ferromagnetic zigzag chains 36] . Since hopping energy is 
minimal for the 180° bond, the competition between the 
hopping and elastic energies generally results in a bond 
angle <p, such that <po < (p < 180° [see Fig. Etc)] . Since <p 
only depends on the nature of the bond (ferromagnetic 
vs. antiferromagnetic), the direction of the oxygen dis- 
placements is the same in all zigzag chains, even though 
neighboring chains have opposite spin. This leads to the 
overall coherent displacement of the center of mass of 
the O atoms with res pec t to the Mn sublattice, similarly 
as proposed in Ref. [la] for the field-induced phase of 
TbMn03. It is easy to see from Fig. [21[b, c) that the 
resulting FE polarization points along the diagonal con- 
necting the next-nearest-neighboring Mn atoms, i. e. the 
a-axis in the orthorhombic setting, in perfect agreement 
with Landau theory. 

The value of P. FigureE^d) shows the calculated abso- 
lute value of the polarization P C aic, defined as the oxygen 
displacement per one "unit cell" containing one Mn and 
two O atoms averaged over MC steps, plotted versus ipo 
for different values of Jaf- P vanishes for ipo = 180° 
for all values of Jaf- In this case both the hopping and 
elastic energies are optimal when (p = (po = 180°, the 
same for ferromagnetic and antiferromagnetic bonds. To 
demonstrate that this is also consistent with the symme- 
try arguments of Landau theory, we note that the sym- 
metry of the lattice with 180° bonds is higher than that 
with tpo < 180°. Particularly, additional inversion cen- 
ters are located at the O sites if (fo = 180°. While the 
E-type magnetic structure breaks inversion symmetry for 
centers located at the Mn sites, it is invariant with re- 
spect to the O-site centers, and thus the magnetic phase 
transition cannot induce FE order. In terms of the Lan- 
dau potential the coefficients c and d must be zero 
in this case. 

The MC calculated value of the polarization P ca i c for 
ipo w 145° corresponding to HoMn03 reaches 0.08. It 
can be shown from the model J2J|, that the FE dis- 
placements scale as i/repE^o- In the simulations we set 
t = 1 and oo = 1, and the polarization in the phys- 
ical units, in the point charge approximation, is given 
by P — ^LoVb ^calo where e is the elementary charge, 
Kp E is the FE phonon stiffnesses in the physical units, 
and Vo is the unit cell volume. Using oq = 2.35 A, Vq = 
226 A 3 0, t = 0.1-0.5 eV k' fe = 1-5 eV/A 2 O, 
we obtain a range of values for P between 0.5 and 12 
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/xC/cm 2 , substantially larger than the P observed in heli- 
cal IMF's. We can also estimate P from the known exper- 
imental results for TbMnOa. Its crystal structure is mod- 
ulated in the collinear sine-wave phase with wave vector 
= 2&m where A;m = 0.29 is the magnetic wave vector. 
The displacement of the O atom n in a Mn-0 chain run- 
ning along the x- or y-direction [see Fig. HJb)] is given 
by 5r n = <5r (— l) n {cos7rfcM + cos[(2n + l)nkM + 2a]}, 
for the magnetic structure defined by S a = S c = 0, 5& = 
S%ax(nk M + a) 0. If k M ? 1/2, then 6r* = 
and the state is paraelectric. The E-phase corresponds 
to fcjvf — 1/2, a — 37r/4, and hence 5r n = 5r repre- 
sents a coherent displacement within one chain. While 
the b and c components of dr^ have different signs in 
different chains, (5rg's are the same. Considering the 
TbMnC-3 value Sr = 1CT 2 A ^ as an estimate for <5rg 
in HoMn03 , we obtain P sw 2 fiC / cm 2 consistent with 
the calculated values. 

Summary. We argue that the symmetry of the spin 
zigzag chain magnetic E-phase in orthorhombic per- 
ovskites with buckling distortion of the oxygen octahcdra 
allows the formation of a polar axis along the a-axis. The 
microscopic mechanism of ferroelectricity is independent 
of spin-orbit coupling, and P can potentially be up to two 
orders of magnitude higher than that in the helical IMF's. 
Note that the ideas described here are general, and the 
E-phase provides just one example. In fact, the mecha- 
nism is similar to the recently proposed explanation for 
the ferroelectricity found in TbMnC>3 at high fields [l5j| . 
and our effort provides a firm theoretical basis for that 
interesting scenario. As to the E-phase nickelates, the mi- 
croscopic model considered here is not applied to them 
directly (3^| . However, we expect that the same interplay 
between the optimization of the electron hopping and the 
elastic energy in the presence of the buckling distortion 
will lead to ferroelectricity. 
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also supported in part by NSF DMR-0443144 and NSF 
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Note added. Following our prediction, B. Lorenz et 
al. |39| reported on experimental evidence of ferroelec- 
tricity in perovskite HoMn03 and YMn03. For the lat- 
ter, the E-phase was proposed in Ref. |l7j . 



[1] T. Kimura et al, Nature 426, 55 (2003). 



[2] N. Hur et al, Nature 429, 392 (2004). 

[3] T. Goto et al, Phys. Rev. Lett. 92, 257201 (2004). 

[4] T. Kimura, G. Lawes, and A. P. Ramirez, Phys. Rev. Lett. 

94, 137201 (2005). 
[5] T. Kimura et al, Phys. Rev. B 71, 224425 (2005). 
[6] G. Lawes et al, Phys. Rev. Lett. 95, 087205 (2005). 
[7] L.C. Chapon et al, Phys. Rev. Lett. 96, 097601 (2006). 
[8] T. Kimura, J.C. Lashley, and A. P. Ramirez, Phys. Rev. 

B 73, 220401 (R) (2006). 
[9] Y. Yamasaki et al, Phys. Rev. Lett. 96, 207204 (2006). 
[10] G.R. Blake et al, Phys. Rev. B 71, 214402 (2005). 
[11] M. Kenzelmann, A.B. Harris, S. Jonas, et al, Phys. Rev. 

Lett. 95, 087206 (2005). 
[12] T. Arima et al, Phys. Rev. Lett. 96, 097202 (2006). 
[13] K. Taniguchi et al, Phys. Rev. Lett. 97, 097203 (2006). 
[14] L.C. Chapon et al, Phys. Rev. Lett. 93, 177402 (2004). 
[15] N. Aliouane, D.N. Argyriou, J. Strempfer, et al, Phys. 

Rev. B 73, 020102(R) (2006). 
[16] A. Mufioz, M.T. Casais, J. A. Alonso, et al, Inorg. Chem. 

40, 1020 (2001). 
[17] J.-S. Zhou and J. B. Goodenough, Phys. Rev. Lett. 96, 
247202 (2006). 

[18] J.L. Garcia-Munoz, J. Rodriguez-Carvajal, and P. La- 

corre, Phys. Rev. B 50, 978 (1994). 
[19] J. A. Alonso, J.L. Garcia-Munoz, M.T. Fernandez-Diaz, 

et al, Phys. Rev. Lett. 82, 3871 (1999). 
[20] M. T. Fernandez-Di'az, J.A. Alonso, M. J. Martinez- Lope, 

et al, Phys. Rev. B 64, 144417 (2001). 
[21] J.-S. Zhou, J. B. Goodenough, and B. Dabrowski, Phys. 

Rev. Lett. 95, 127204 (2005). 
[22] H. Katsura, N. Nagaosa, and A.V. Balatsky, Phys. Rev. 

Lett. 95, 057205 (2005). 
[23] LA. Sergienko and E. Dagotto, Phys. Rev.B 73, 094434 

(2006). 

[24] M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). 

[25] M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005). 

[26] C. Ederer and N.A. Spaldin, Curr. Opin. Solid State 
Mater. Sci. 9, 128 (2005). 

[27] B. Lorenz, Y.Q. Wang, Y.Y. Sun, and C.W. Chu, Phys. 
Rev. B 70, 212412 (2004). 

[28] S. Goshen, D. Mukamel, H. Shaked, and S. Shtrikman, 
Phys. Rev. B 2, 4679 (1970). 

[29] A. B. Harris, T. Yildirim, A. Aharony, and O. Entin- 
Wohlman, Phys. Rev. B 73, 184433 (2006). 

[30] Although the E-phase is a special case of the simple 
sine- wave structure, the continuum theory considered in 
Ref. [2^1 . which predicts P = 0, is not applicable to mod- 
ulation vectors at the edge of the Brillouin zone. 

[31] E. Dagotto, T. Hotta, and A. Moreo, Phys. Rep. 344, 1 
(2001). 

[32] T. Hotta et al, Phys. Rev. Lett. 90, 247203 (2003). 

[33] D.V. Efremov, J. van den Brink, and D.I. Khomskii, Na- 
ture Materials 3, 853 (2004). 

[34] J.C. Slater and G.F. Koster, Phys. Rev. 94, 1498 (1954). 

[35] T. Kimura, S. Ishihara, H. Shintani, et al, Phys. Rev. B 
68, 060403(R) (2003). 

[36] Although virtual hopping is allowed along the zigzag 
chains, the E-phase is a band insulator [32^|. 

[37] T. Kimura, private communication. 

[38] T. Hotta and E. Dagotto, Phys. Rev. Lett. 92, 227201 
(2004). 

[39] B. Lorenz, Y .Q. Wang, and C.W. Chu, 
cond-mat/0608195 



